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FURTHER  DEVELOPMENT  OF  HIERARCHICAL 
PLASTICITY  MODEL  FOR  DISCONTINUITIES 


1 


This  model  represents  a  generalized  and  modified  development  of  the 
model  presented  in  Part  A  of  this  report.  The  proposed  generalized  model 
allows  for  detailed  consideration  of  the  normal  and  shear  components  of 
the  joint  response,  and  is  verified  for  both  normal  stress  and  normal 
stiffness  controlled  conditions. 

STRESSES  AND  DISPLACEMENTS 
Stresses 

A  joint  is  usually  idealized  as  a  planar  thin  layer  and  only  the 
stresses  on  this  plane  are  considered  (Fig.  1).  Thus,  only  the  normal 
stress  and  the  shear  stress  t  are  included  in  the  modeling  (Fig.  2)  , 
and  the  influences  of  the  other  stress  components  are  ignored.  Since 
due  to  the  roughness  of  the  joint,  there  are  only  finite  number  of  con¬ 
tacts  and  stresses  are  concentrated  at  these  contacts,  an  assumption  is 
also  made  that  the  stresses  on  the  joint  plane  are  uniformly  distributed 
so  the  stresses  are  nominal  rather  than  actual. 


Normal  Displacement 

Under  normal  compression,  the  joint  will  be  compressed,  and  under 
tension,  the  joint  can  separate.  During  shear,  dilation  may  occur 
(Fig.  3),  and  this  dilation  is  caused  by  the  sliding  of  the  asperities 
on  the  opposite  walls  of  the  joint.  These  compressions,  separations 
and  dilations  are  called  normal  displacements  since  they  occur  in  the 
normal  direction  to  the  joint  plane. 


The  total  normal  displacement  v  can  be  decomposed  into  three  parts: 


the  elastic  and  plastic  deformation  of  the  asperities  v  and  v  ,  re- 

e  p 

spectively,  and  the  slip  displacement  v  ,  as 


V 


=  V  +  V  +  V 

e  p  s 

Cl) 

=  V  +  V 

e  r 


where  is  called  the  relative  normal  displacement. 

The  normal  slip  displacement  v  is  the  projection  of  the  sliding 
and  slip  displacement  on  the  normal  direction.  In  a  normal  compression 
test  of  an  initially  mated  joint,  the  normal  compression  is  mainly  due 
to  the  deformation  of  the  asperities  and  the  normal  slip  displacement  is 
negligible.  On  the  other  hand,  however,  in  a  shear  test  with  constant 
normal  stress,  the  vertical  slip  displacement  is  a  major  part  of  the  di¬ 
lation  and  the  elastic  and  plastic  deformations  of  the  asperities  are 
relatively  small  compared  to  the  slip  displacement. 


Similar  to  the  normal  displacement,  the  tangential  displacement  u 
can  be  decomposed  into  three  parts: 


u  =  u  +  u  +  u  (2) 

where  ug,  represent  the  elastic  and  plastic  shearing  deformations  of 
the  asperities,  respectively,  and  is  the  tangential  slip  displacement. 

It  is  very  difficult  to  separate  the  plastic  deformation  u^  and  the 
slip  displacement  u  ,  so  they  can  be  combined  as  relative  tangential  dis¬ 
placement  u  ,  as 
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Usually,  the  plastic  deformation  is  quite  small  compared  to  the 
slip  displacement  u^ .  In  some  cases,  the  plastic  deformation  of  the 
asperities  u^  is  ignored  in  modeling. 

Plasticity  for  Joints 

The  plasticity  theory  was  developed  for  continua.  In  the  plasticity 
theory,  there  is  a  yield  function,  F,  that  serves  as  a  field  equation, 
and  a  flow  rule  which  determines  the  plastic  deformation  directions. 

For  joints,  however,  plasticity  can  not  be  directly  used.  Firstly, 
the  joint  is  two-dimensional ;  hence,  a  special  form  of  the  yield  function 
is  needed.  Secondly,  the  irrecoverable  displacements  are  not  necessarily 
plastic  deformations,  but  a  combination  of  sliding,  slipping  and  plastic 
deformations  of  the  asperities.  Thus,  the  displacement  has  a  preferred 
direction  which  is  determined  by  the  surface  geometry  of  the  joint.  The 
plastic  deformations  of  the  asperities  enter  as  a  "correction"  to  this 
preferred  direction.  When  the  joint  is  relatively  hard  or  the  stress 
level  is  relatively  low  such  that  plastic  deformations  are  very  small 
compared  to  the  slip  displacement,  we  can  neglect  such  "correction." 

By  specializing  the  general  three-dimensional  yield  function  for 
continua  developed  by  Desai  et  al.  [1],  a  yield  function  F  for  joints  is 
obtained  by  Desai  and  Fishman  [2],  Fishman  and  Desai  [3]  and  Fishman  [4] 


^  2  n  2 

F=x  +  a  a  -  y  a 

n  n 


where  n,  y  are  response  functions  and  a  is  the  hardening  or  growth  func¬ 


tion.  The  hardening  function  a  can  be  expressed  as 
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where  a^,  are  the  normal  and  shear  hardening  functions,  respectively. 

The  flow  rule  for  joint  is  different  from  that  of  continua  due  to 
the  existence  of  the  slip  displacement.  The  flow  rule  for  the  plastic 
deformation  for  the  two-dimensional  case  can  be  expressed  as 


=  X 


3a 


where  d  denotes  increment,  Q  is  the  plastic  potential  and  \  is  the  scalar 
proportionality  parameter. 

For  the  joint,  the  plastic  deformation  is  the  difference  between  the 
relative  displacement  and  the  slip  displacement,  Eq.  (1).  Thus, 


dv  =  dv  +  dv 
r  s  p 


du  =  du  +  du 
r  s  p 


dv 

s 

du 

s 


+  \ 


3Q_ 

3a 


3Q 

3x 


(7; 


where  dvr,  duf  are  the  incremental  normal  and  the  tangential  relative 
displacements,  respectively,  and  dvs,  dug  are  the  incremental  normal  and 
tangential  slip  displacements,  respectively. 

If  the  plastic  deformation  is  negligible,  then  the  relative  dis¬ 


placement  is  approximately  equal  to  the  slip  displacement. 


Stress-Displacement  Relations 


Normal  Compression  Test 


In  the  normal  compression  test,  for  initially  mated  joints,  the 

shear  stress  x  and  the  tangential  displacement  are  equal  to  zero.  Thus, 

the  yield  function  F  reduces  to  F  : 

n 

F=aan-ya2=0  (I 

n  n  n  n  '■ 


and  from  the  consistency  condition,  dF  =  0; 


3F  3F 

dF  =  -*■ — —  do  +  ^  dv  s  0 
n  3a  n  3v  p 
n  P 


there,  the  normal  displacement  is  mainly  deformation  of  the  asperities 


and  the  slip  displacement  is  assumed  to  be  zero.  So, 


dv  =  dv  +  dv 
e  p 


The  elastic  vertical  displacement  dvg  is  given  by 


dv  =  — 
e  K 


where  K  is  the  unloading  elastic  normal  stiffness  of  the  joint.  The 
n  &  J 

plastic  normal  displacement  is  obtained  from  Eq.  (9)  as 


3F  3F 

dv  =  -  [  0—^)  /  (V-^)  ]  do 
p  L  v3a  ^3v  J J  n 

r  n  p 


So  the  total  displacement  is 
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(15) 


dv  =  [(1T  '  dan 

n  n  p 

Different  models  can  be  obtained  by  choosing  different  hardening 
function  a^.  One  possible  model,  based  on  observed  laboratory  behavior, 
can  be  expressed  as 


an  =  y  Exp  [  (2-n)  (hx  vp  +  h^]  (14) 

where  y,  n,  h1  and  h_,  are  parameters  and  v  is  the  normal  plastic  compression. 
1  4  p 

Shear  Tests 

In  the  previous  section,  the  normal  compression  behavior  (or  the  in¬ 
cremental  normal  stress  and  normal  compression)  is  modeled  by  using  the 
yield  function  F  .  During  shear,  if  we  assume  that  the  normal  stress  and 
normal  plastic  deformation  relation  is  not  changed,  then  we  can  assume  that 


n  2 

F  =  a  a  -ra  =0 
n  n  n  n 


(15) 


still  holds.  The  yield  function  relevant  to  the  shear  condition  can  be 
expressed  as,  Eq.  (4), 


F 

T 


+  a 


a  a 

n  n 


-  y  a 


By  applying  Eq.  (15) ,  we  obtain 

F  =t2+y(cc  -  1)  a  2  =  0 
t  x  n 


(16) 


(17) 


where  a  is  the  hardening  or  growth  function  and  can  be  expressed  as 


1  n4 

=  1  -  -  Exp  (h,)  x  ur 


(18) 
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h,  =  u  a  +  b 
3  n 


h  =  c  a  +  d 
4  n 


where  a,  b,  c,  d  and  y  are  constants,  and  is  the  relative  tangential 
displacement . 

The  total  incremental  displacements  can  now  be  expressed  as 


du  =  du  +  du  =  du  +  — 
r  e  r  K 

T 


dv  =  dv  +  dv  =  dv  +  — 
r  e  r  K 


where  ,  Kn  are  the  unloading  shear  and  normal  stiffnesses,  respectively. 

Because  the  plastic  deformations  of  the  asperities  are  in  many  cases 
very  small  compared  to  the  slip  displacement.  The  flow  rule  can  be  approxi¬ 
mately  assumed  as 


dv  =  s  du 
r  r 


where  g  is  the  geometrical  function,  as 


where 


f  =  [1  -  Exp  (-k  u^.  )  ] 


where  v  is  the  ultimate  dilation  and  k,  m  are  constants, 
u 


% 


Shear  tests  can  be  classified  as  (a)  constant  normal  force  stress 
and  (b)  constant  normal  stiffness.  However,  very  often  both  situations 


may  be  coupled. 

For  a  shear  test  with  a  constant  normal  loading  stiffness,  category  (b) 
k^,  (or  a  CNK  test),  when  there  is  dilation  dv,  the  loading  system  will 
increase  the  normal  stress  such  that 

da  =  K,  dv  (25) 

n  d 

This  test  is  denoted  by  CNK  here.  Shear  test  with  constant  normal  stress, 
category  (a),  as  denoted  by  CNS. 

From  the  consistency  condition: 

3F  3F  3Ft  3F 

dF  =  dx  +  -7T- —  da  +•  ■  du  +  -t— -  dv  =0  (26) 

T  9t  3a  n  3u  r  dv  r  v 


and 


da  = 


Kd  g 


n  K 


K .  ^du  "  K  ^ 
d  T 


(30) 


PARAMETERS 

The  parameters  to  be  determined  are  the  elastic  stiffness,  the  con¬ 
stants  in  the  hardening  functions,  the  constants  in  the  geometrical 
function  and  the  constant  for  the  ultimate  state. 

Elastic  Constants 

The  normal  compression  elastic  stiffness  and  the  shear  elastic 
stiffness  K  are  the  elastic  constants  used  in  the  model.  They  are 
assumed  equal  to  the  unloading  stiffness  in  the  normal  compression  test 
(Fig.  4)  and  the  shear  test  (Fig.  5)  ,  respectively. 

Hardening  Functions 

The  constants  in  the  normal  hardening  function  a  are  h,  and  h_. 

n  1  2 

They  are  obtained  from  the  in  plot  (Fig.  6),  and 

In  a  =  h.  v  +  h_  (31) 

n  1  p  2 

The  parameters  in  the  shear  hardening  function  a  are  h,  and  h  .  They 

T  J  4 

are  obtained  from  the  in  (— )  -  in  plot  (Fig.  7)  and 

n 

in  (— )  =  h7  In  u  +  h.  (52) 

'■a  3  r  4 

n 

and  h,,  h^  are  further  determined  by 
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Normal  Stress 


Shear  Stress 


Fig.  5  Unloading  Elastic  Stiffness  K  from  Shear  Test 
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Normal  Plastic  Displacement 


Fig.  6  Plot  to  Find  Constants  and 


« 
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Fig. 


£n 


Plot  to  Find  Parameters  h  and  h 
Which  Leads  to  Find  Constants  a,  b. 
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[ 


h  ,  =  a  o  +  b 

3  n 

(33) 

h.  =  c  a  +  d 

4  n 


where  a,  b,  c,  d  are  constants. 


Geometrical  Function 

To  determine  the  geometrical  function,  three  constants  are  needed: 
v  ,  k  and  m,  and  the  geometrical  function  g  is  expressed  as 


g 


(34) 


where 


% 


f  ■  vu  [1  -  Exp  (-k  u/)] 


(Fig. 


v 

These  constants  are  obtained  from  the  Jin  [Jin  ( - - - )] 

u  r 

8),  where  v^  is  the  relative  normal  displacement. 


(35) 

&n  u^  plot 


Constant  of  Ultimate  State 

At  the  ultimate  state,  a  =  0,  so  the  yield  function  becomes 


F 


f 


=  0 


(36) 


so 


2 


a 

n 


and  y  is  obtained  by  best  curve  fitting  to  the  Df  ultimate  stress  states 
for  different  normal  stresses  (Fig.  9a). 
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£n  [£n 


Shear  Stress 


V 


m 


% 


% 


% 


Constant  n 

The  constant  n  is  usually  determined  at  the  phase  change  point. 

Fig.  9(b),  and  for  a  joint,  this  point  is  often  hard  to  select,  .as  this 
time,  based  on  previous  tests  [4],  n  =  2.5  is  adopted.  However,  later 
it  would  be  further  investigated. 

BACK  PREDICTIONS  AND  VERIFICATION 

In  order  to  verify  the  model,  three  sets  of  test  data  are  used.  The 
first  data  set  is  from  a  paper  by  Yoshinaka  and  Yamabe  [5].  The  joint  is 
artificially  made  by  cutting  an  originally  intact  rock  with  the  diamond 
saw.  The  tests  performed  are  the  normal  compression  test  and  the  shear 
test  with  constant  normal  stresses. 

The  second  data  set  is  from  the  tests  performed  by  Fishman  [4] .  The 
joint  is  prepared  by  casting  the  concrete  against  a  saw-shaped  metal  sheet. 
Here,  only  the  shear  tests  with  constant  normal  stresses  are  used  in  back 
predictions . 

The  third  data  set  is  from  Dight  and  Chiu  [6].  The  joints  are  very 
rough  surfaces  of  rock  on  concrete.  These  shear  tests  involve  constant 
normal  stiffness. 

The  above  three  tests  are  designated  as  normal  compression  test  (NCT) , 
the  shear  test  with  constant  normal  stresses  (CNS)  ,  and  the  shear  test 
with  constant  normal  stiffness  (CNK) . 

For  back  predicting  the  NCT  tests,  the  constants  required  are  Kn, 

h^  and  h-,.  For  CNS  and  CNK  tests,  the  constants  required  are  Kx ,  y,  a,  b, 

c,  d  and  v  ,  k,  m. 
u 
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Table  1  gives  the  constants  for  all  three  data  sets.  Table  2  indi¬ 
cates  all  the  figures  (from  Fig.  10  to  Fig.  17)  that  show  the  correlation 
between  the  data  and  the  constants  found.  For  all  the  data  sets,  the  corre¬ 
lation  is  satisfactory. 

Back  predictions  are  made  and  compared  to  the  data  for  all  three 
data  sets  (Fig.  18  -  Fig.  23). 

Figures  18  and  19  back  predict  the  relations  of  normal  stress  and 
normal  displacement  v  and  the  relation  of  shear  stress  x  and  tangential 
displacement  u  for  data  set  #1.  The  results  are  very  satisfactory. 

Figures  20  and  21  give  the  back  prediction  of  the  relation  of  shear 
stress  x  and  tangential  displacement  u  and  the  relation  of  dilation  v  and 
tangential  displacement  u  for  data  set  #2.  This  is  a  shear  test  with 
constant  normal  stresses.  But  under  different  normal  stresses,  the  dila¬ 
tions  are  almost  the  same  since  the  joint  is  very  stiff  and  the  dilations 
are  mostly  rigid  slips. 

Figures  22  and  23  give  the  back  predictions  for  data  set  #3.  Because 
this  is  a  shear  test  with  constant  normal  stiffness  CNK,  the  stress- 
displacement  relations  are  quite  different  from  the  CNS  tests.  A  very 
interesting  point  is  that  while  the  normal  stiffness  keeps  constant,  the 
normal  stress  a  increases  and  the  shear  stress  x  decreases  after  it 
reaches  a  peak.  This  is  a  complicated  situation  since  the  softening  occurs 
in  shear  but  strengthening  or  hardening  happens  in  normal  direction.  The 
model  catches  this  feature  satisfactorily. 

Summary  and  Conclusions 

A  plasticity  model  for  joints  and  discontinuities  is  presented.  Due 
to  many  special  features  of  the  joints  and  discontinuities,  the  plasticity 
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Fig.  18  Back  Prediction  -  Normal  Stress  o 
and  Normal  Displacement  v  Relation 
in  Normal  Compression  Test  (Data  Set  #1) 
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Tangential  Displacement  u  (mm) 


Fig.  19  Comparisons  Between  Back  Predictions  and 
Observations  for  Shear  Stress  and  Shear 
Displacement  Relations  Under  Various  Normal 
Stresses  a  in  Shear  Tests  (CNS)  (Data  Set  #1) 


Observed 


Fig.  20  Comparisons  Between  Back  Predictions  and  Observations 
for  Shear  Stress  and  Tangential  Displacement  Relation 
(Jnder  Various  Normal  Stresses  a  in  Shear  Test  (CNS) 
(Data  Set  #2)  n 


Normal  Displacement  v  (inj 


Fig.  21  Comparisons  Between  Back.  Predictions  and  Observations 
for  Dilation  and  Tangential  Displacement  Relation  in 
Shear  Test  (CNS)  (Data  Set  #2) 
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23  Comparisons  Between  Back  Predictions  and  Observations 
for  Dilation  and  Tangential  Displacement  Relation  in 
Shear  Test  (CNK)  (Data  Set  #3) 
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theory  is  applied  with  modification.  The  major  modification  is  made  for 
the  flow  rule  since  the  direction  of  displacement  is  determined  mainly 
by  the  geometry  of  the  joint,  especially  for  very  stiff  joints. 

Verification  is  made  for  three  typical  tests,  the  normal  com¬ 
pression  test  (NCT) ,  the  shear  test  with  constant  normal  stress  (CNS) 
and  the  shear  test  with  constant  normal  stiffness  (CNK) .  The  back  pre¬ 
dictions  are  very  satisfactory. 
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